Surface tension of liquids near the critical point 
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A simple analytical microscopic expression for the surface tension of liquids 7 is obtained which 
is in a good agreement with available data of numerical experiments. We apply the integral trans- 
formation that maps the fluid Hamiltonian onto the field-theoretical Hamiltonian and show that 
the order parameter of the effective Hamiltonian corresponds to the one-body potential in the fluid. 
Revealing the physical meaning of the order parameter allows calculation of 7. 



The modern theory of nonuniform fluids is based 
mainly on the integral equations hierarchy or density 
functional theory [Q. Recent applications of the latter 
theory to the problems of interface formation and sur- 
face tension demonstrated its accuracy and effectiveness 

H . Although the density functional theory often uses 
ad hoc approximations, rather than regular expansion 
schemes |^ , it is essentially the first-principle theory: No 
model parameter are involved, and for a given interparti- 
cle potential, temperature, etc., the density distribution 
and surface tension may be obtained ^ . Being phys- 
ically very clear and straightforward, the density func- 
tional theory is not very easy in practical applications, 
since sophisticated numerical analysis is required ^. 
It does not generally provide any analytical expression for 
the surface tension. The other limitation of the density 
functional theory refers to its strictly mean-field nature 
Apart from its phenomenological version |^ (where 
model parameters were used) the density functional the- 
ory is not suited to study the critical region, where the 
critical fluctuations play a clue role. 

In somewhat alternative, field-theoretical approach 
an order parameter (f>{r) is introduced and functional de- 
pendence of the free energy (or of the effective Hamilto- 
nian) on (j){r) is postulated. For the most frequently used 
square-gradient model it reads (e.g. |||, ^) 

Pn = I dr[(«;/2)(V0)2 + T/(0)] (1) 

where the first term accounts for the free energy penalty 
due to inhomogeneity and V{(j)) is chosen to mimic 
a possibility of the two-phase coexistence. Its sim- 
plest choice corresponds to the Landau-Ginzburg- Wilson 
(LGW) Hamiltonian: 

y(0) = {a2/2l)cl>\r) + M4l)cb\r) - h{v)(p{v) , (2) 

with the free energy, —j3F = log/ l?0exp(— /3H) (/3 = 
l/fc^T), which accounts for all possible distributions of 
the order parameter, and not only for the extremal one, 
as in the density functional theory. The powerful field- 
theoretical methods, elaborated especially for the LGW 
Hamiltonian |6j may be then applied to describe the 



critical fluctuations, and a simple analytical expression 
for the surface tension may be derived |^ . 

The order parameter 0(r) in (|l|,^ may be of any phys- 
ical nature, (density, composition of the fluid mixture 
0, magnetization etc.), while the coefficients of the 
Hamiltonian Ti are some model parameters for which mi- 
croscopic expressions are not provided. Therefore the 
field theoretical approach is not now, strictly speaking, a 
microscopic theory. 

The present study is addressed to overcome this flaw. 
We obtain an effective field theoretical Hamiltonian for 
the fiuid with microscopic expressions for its coefficients. 
We show that the order parameter of this Hamiltonian 
corresponds to the one-body microscopic potential, and 
its average to the first-order direct correlation function. 
Revealing the physical nature of the order parameter we 
obtain an analytical expression for the surface tension 
which is in a good agreement with numerical experi- 
ments. 

Referring for detail to ||] we briefly sketch derivation 
of the effective LGW Hamiltonian, which is similar to 
that of Hubbard and Schofield . We start from the 
fiuid Hamiltonian H = Hn + Ha + H^^: 

i<j i<j j 

where Vr{r) denotes the repulsive part of the interac- 
tion potential, —v{r) denotes the attractive part and 
g{r) denotes the external potential; {rj} are coordi- 
nates of the particles and = — rj. The last 
two terms of the fluid Hamiltonian (^ may be writ- 
ten using the Fourier transforms of the density fluctu- 
ations, rik = (l/v^) X^^i e"*'' ''' , of the attractive po- 
tential, Vk = J v{r)e~'^^''^dr and of the external potential 
9k = J g{r)e^^^''^dr as 

- (1/2) ^'fcrikn-k + (l/2)i;(0)iV + 5kn_k , (4) 

k k 

where N is the number of particles, fl = L'^ is the volume 
of the system {L — > 00), and summation over ki — 
with I = x,y,z, and ni ~ 0,±1,... is implied. Let /i 



2 



be the chemical potential of the system with the "com- 
plete" Hamiltonian, and be the chemical poten- 
tial of the "reference" system, with the Hamiltonian, Hji, 
which has only repulsive interactions. Then the grand 
partition function, S, may be expressed in terms of the 
grand partition function, of the reference fluid as Q 



^ — 



exp < + /3 ^ [ Y?^k?^-k - "kg- 




where /i' = /i — /i + ^^(O) and ( )^ denotes an average 
over the reference system at temperature T with chemical 



potential iir. Using the identity: exp(ia 



2^2 



bx) 



(27ra2)-i/2 J+^exp[-i(y + bf/a^ + xy]dy, we obtain 
for the ratio Q = S/S^;: 

Q oc y"]^d(/)k ^exp |^0k»^-k|^ cxp|^fi^/^0o 

- ^ E ^fc ' (</>k + P9k) (<^-k + /35-k) I ; (6) 

here integration is to be performed under the constrain 
= (p'^, and a factor which does not affect the subse- 
quent analysis is omitted. Applying the cumulant theo- 
rem to the factor (exp {X^k 0k'T--k})^ one can write ||ic[| 



exp{—f3H) 



(7) 



"=2 ki,...k„ 

where the coefficients of the effective field theoretical 
Hamiltonian Ti read for glr) = 0: h = fJ-'v^^ + p, and 

{t2(ki,k2) {l/2\)Sk,+k^^o{P^^v^^^ - (nkin-ki)ci?j , 
M„ (ki, . . . k„) = -(r!"/2-i/^,) ^^^^ . . . ^^j^^ 

for n > 3. Here ( denotes the cumulant average 
calculated in the (homogeneous) reference system and 
p = < no >cR^ N/Q is the fluid density. Accord- 

ing to (^, Q has the form of a partition function. 

As it follows from (||) the coefficients of Ti. depend on 
correlation function of the reference fluid having only 
repulsive interactions. Using definitions of the parti- 
cle correlation functions of fluids one can express 
the cumulant averages (nk^ • • • "k„)cfli ^^id thus the co- 
efficients {in (ki, . . . k„) in terms of the Fourier trans- 
forms of the connected correlation functions hi,h2, . ■ .hn 
of the reference system. These are defined as /ii(ri) = 
(5(ri), /i2(ri,r2) = g2(ri,r2) - 1, /i3(ri, r2, ra) = 
53 (r 1 , r2 , rg ) - .92 (r 1 , r2 ) - 32 (r 1 , rs ) - 32 (r2 , rs ) + 2 , etc. , 



where gi{ri,...ri) are Z-particle correlation functions 
In particular, the first few coefficients read M 



W2 = (l/2!) {f]vo)-'-p-p^h2{-ki) 



Jki+k2,0 



(9) 



U3 



-(1/3!) { 



P + P 



/l2(kl) +/l2(k2) +/l2(k3) 



+p3/i3(ki;k2)}(5ki 



+k2+k3,0 , 



(10) 



where hi are the Fourier transforms of hi. 

To obtain the conventional square gradient form of the 
effective Hamiltonian we perform the small-fc expansion 
of the coefficients Un (ki, . . .k„). From the structure of 
the LGW Hamiltonian (|l|), which has the only gradi- 
ent term (V0)^ ^ fc^'/>k<?!>-k, follows that only {t2 should 
be expanded as U2 — "2(0) — W2(0)"fc^ + • • while the 
other coefficients iin, n > 3 should be taken at zero wave- 
vectors, as {t„ (0, 0, . . . 0) 1^. Thus, as it is seen from 
(|,0) only /i2(0), hiOy and hi{0) = hi{0,0, . . .0), I > 2, 
are needed. To find the latter we use the chain relation 
for the successive correlation functions H 



g 

zop-gj^P^hi{0) 



lhiiO) + hi+iiO) 



(11) 



expressing each /i/+i(0) in terms of /i;(0) and its den- 
sity derivative. Here zq = f3~^ (dP/dp)^^ is the reduced 
isothermal compressibility and P is the pressure. With 
hi{0) = 1, Eq.jTll) allows to express iteratively each of 
hi(0), (which refer to the reference fluid) in terms of the 
reduced compressibility of the reference fluid zq and its 
density derivatives dz^/dp, d'^zo/dp^, etc. 

For the reference system with the only repulsive in- 
teractions one can use the hard-sphere system with an 
appropriately chosen diameter [ pT| . For soft (not im- 
pulsive) repulsive forces a simple relation [ pT| d = 
/J' [1 — exp (— /3u,.(r))] gives the effective diameter of the 
hard-sphere system, corresponding to a repulsive po- 
tential Vr{r) vanishing at r > a. The fairly accurate 
Carnahan-Starling equation of state for this system 
yields for the reduced compressibility 



zo = (1 - ?])** / {1 + 4t] + V - 4??^ + V^) 



(12) 



where 77 = ■^d'^p. For the hard-sphere reference sys- 
tem one can also find /i2(0)". This may be done ex- 
pressing /i2(fc) in terms of the direct correlation func- 
tion C2(fc), as /i2(fc) — C2{k) / [1 — pc2{k)] expand- 
ing C2(fc) as C2(fc) = £2(0) — C2(0)"fc^ + • • and using 
52(0)" = -(7rdVl20)(l6-ll7?-h 4772)(l-?7)"'' which 
may be obtained from the Wertheim-Thiele solution 
for the direct correlation function ||] . 

Following these lines we find explicite expressions for 
all coefficients of the effective Hamiltonian (|^). To recast 
this into conventional form we (i) perform trans- 

formation from variables 0k to the space-dependent field 
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(/)(r), (ii) omit all terms in (|^) with powers of the field 
higher that fourth, and (iii)remove the cubic term with 
respect to the field, which V^cj)) does not contain. This 
may be done by the shift — > cf) + (j), with chosen to 
make the cubic term vanish. As the result, we arrive 
at the effective Hamiltonian (|l|) with microscopic expres- 
sions for its coefficients. The coefficients of V{(p) read 

U4 = -pzo[zf + zo{zo + 4zi + Z2)] 

a2 ^ ilSvo)-^ ~ p{zo + zy2ui) (f3) 

h = p'vq^ + p{a2 + zl/6u4)z3/u4 + p . 

with zq given by (p^), zi = pdz^/dp, z^ = p^ d'^ zq / d p^ 
and Z3 = ~pzo{zo + zi), while the coefficient k reads 



(3/407rd) [Xls/P^off - B] 



(14) 



where B = 4r/2(l - r/)**(f6 - Ut] + 4ry2)/(l + Arj + Ajf ~ 
Arf + rf)"^ and constants Ccff and Acff (which appear due 
to the expansion of U2) characterize the effective depth 
and effective width of the attractive part v{r): e^s — 
{4T:d^ J v{r)dr and A^^j = (Svod"^ J v{rydr. 

To this point, we have derived the effective field theo- 
retical Hamiltonian with explicit microscopic expressions 
for its coefficients. Now we show that the order parame- 
ter 0(r) which has been introduced as a "technical" vari- 
able of the integral transformation has a clear physical 
meaning: it is related to the local microscopic potential. 

To show this we notice that from (||) it directly follows 
91ogS/9(7„k = 91ogQ/9g_k = — /5(?ik); on the other 
hand Eq.(||) allows to write for 5 — > 0: dlogQ/dg-u = 
~'^k^ (0k) , where the averaging is to be understood as the 
integration over all distributions of the order parameter. 
Recalling that we have shifted the field by the constant 
(f>, we conclude that ((/)k) = I3vk{ni^) + ^<^k,Oj or 



{(f>{r)) = P / p(ri)w(r - ri)dri + 



(15) 



where p(r) = (n(r)) is the average one-particle density. 
Eq.(|l5|) shows that (up to cj)) the average order parame- 
ter at point r equals to the effective one-body potential 
(in units of fcsT), therefore, — ((/)(r)) has a similar phys- 
ical meaning as c^^^(r), the first-order direct correlation 
function in fluid theory One may thus conclude that 
the order parameter (j)(r) corresponds to the microscopic 
one-body potential within the fluid, the effective Hamil- 
tonian gives energy, written in terms of the microscopic 
potential, and the partition function equals to the (func- 
tional) integral over this field (see (0)). It is also worth 
to note that ((/>(r)) changes more smoothly in space than 
p{r), and that for homogeneous fluids {p{r) = p = const) 
((/)(r)) = f3vop + (f> = const. 

Now we calculate the surface tension using the mean 
field approach where only the extremal field 0* (r) , which 



minimizes the free energy, is taken into account, i.e. 
Fmf = Ti-icj}*)- For the flat interface with V = d/dx, 
the equation for the extremal field reads |P, 0] 



Kd^(j)*/dx'^ = dV{(p*)/d(j)* 



(16) 



In the bulk of the two phases, i.e. far from the sur- 
face, the order parameter takes constant values, (j)^ at 



and 



at X 



00, which are related to 



the mean densities of these phases: the liquid, pi, and 
the vapor, pg density. Using the fact that in the bulk 
'f'1.2 — PvoPi,g + one can show that 4>\ and 4>2 may 
be obtained by the standard double-tangent construe- 
tion, V'm = and V ^ + ^^V {4>*i) - ^'(-^2) + 

(t>2V'{4>2)- If ^6 choose interface located at x = 0, ( p^ ) 
yields \K{d(t)*/dxY = V{(j)*) - V{(t>l) for a; < and 
\K{d(j)*/dxf = V{(j)*) - V{(l)*2) for a; > 0. The surface 
tension is equal to the difference per unit area between 
the free energy, calculated for the space-dependent (j)* (r) 
and that for (j)^ for a: < and ^2 for a; > 0. If the or- 
der parameter at the interface equals (p^, which may be 
chosen from the condition </)J < 0q < ^'(<Ao) = the 
surface tension reads (with Vi^2 — y{<f>i,2))'- 



(17) 

Now we choose the system for which the coefficient 
h in ( |l3| ) vanishes, i.e. h = 0. Then for V = 



^a2(, 



iu4(/>*'* we obtain 



±(-6a2/w4)l/^ 



= 0, and the symmetric solution to Eq.(|iq), 0*(x) = 
(-602/^4)'/" tanh(a:/^o) (Co =_ ^-K/2a2)~§, 0, with 
zero over the volume average, 0* = ft^^ J 4i*{r)dr = 0. 

Averaging over the volume (p^), yields (p* — (3pvo + (f>, 
which implies that cj) = —(3voi5, with p — fl"^ J p{r)dr = 
N/D, being the averaged over the volume density. Since 
0* = —02 and simultaneously 0i 2 = '^P'VoPi,g + 4>^ 
finally conclude that p — {pi + Pg)/"^, i-S- that the av- 
eraged density of our system is the mean between the 
liquid and vapor density. Naturally, this is the density of 
our homogeneous reference system, which has the same 
volume and number of particles. With the above 2 
and 03? integration in ([l7|) is easily performed yielding: 



a) 



(18) 



where microscopic expressions for the constants 02, U4, 
are given by Eqs.(^ 14) in which the density p — 



{pi + Pg)/2 of the reference fluid is to be used. 

Not far from the critical point (pc, T^), one can ap- 
proximate, (pi + Pg)/2 ~ pc and thus use pc as the 
reference density. In particular one can write for a2'. 
a2 ^ a2iP, Pc) = iPvo)-^-Pcizo + zl/2ui)c (see (|3[)). If 
we then use the the mean field condition for the critical 
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point, a2{Pc,Pc) = 00, 
tension: 



we obtain a simple expression: 
- ar, and finally for the surface 



(19) 



where a = {(3cVq) ^, t — {Tc — T)lTc., and coefficients 
M4^c and Kc are to be calculated a,t p = pc, T = Tc- 
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FIG. 1: Reduced surface tension 'ja^/ksT as a function of 
the reduced temperature ksT/e for the Lennard- Jones (LJ) 
and hard-core Yukawa (HCY) fluid. Curves - theory, Eq.([f9|): 
1 - LJ-fiuid; 2, (2'), 3, 4 - HCY-fluid for A = 1.8, A = 3.0 and 
A = 4.0 correspondingly. Points - numerical data: diamonds 

- LJ (MD circles - LJ (MC O), stars and triangles up 

- HCY for A = I.SJMC and MD 0), triangles down - HCY 
for A = 3.0 (MD Q), squares - HCY for A = 4.0 (MD 0). 
Critical parameters are taken from [^] for the LJ-fluid, from 
m (curves 2-4) and from ||l|] (curve 2') for the HCY-fluid. 
a, e and cr, e, A are parameters of the LJ and HCY potential 
lHll4|. 



The theoretical expression for the surface tension ( pj| ) 
has been compared with available data of numerical 
experiments for the Lennard-Jones (LJ) and hard-core 
Yukawa (HCY) fluids, for which the standard WCA par- 
tition (see e.g. ||ll|) of the potential into attractive and 
repulsive parts has been applied ||^. As follows from 
Fig. 1 our theory is in a good agreement with the numer- 
ical experiments, except for the very close vicinity of the 
critical point, where the mean field theory loses its accu- 
racy. Eq. (hof) is quite sensitive to the critical parameters 



Tc- While these are known quite accurately for the 
LJ fluid, they are estimated with much larger uncertainty 
for the HCY fluid. This is shown in Fig.l where two the- 
oretical curves (2 and 2') correspond to the same HCY 
fluid but with pc, T^, taken from different references (pc 
and differ by about 4%). 

In conclusion, we derive an analytical microscopic ex- 
pression for the surface tension which allows (for the 
first time to our knowledge) to obtain this quantity for 
the given temperature, density and interaction potential. 
This relation is in a good agreement with numerical ex- 
periments in the near-critical region. Our theory is based 
on the effective field theoretical Hamiltonian, which or- 
der parameter corresponds to the microscopic one-body 
potential and its average to the first-order direct correla- 
tion function. Although our mean field result is not valid 
in the very close vicinity of the critical point, its gener- 
alization to account the critical fluctuations (e.g. within 
the one-loop approximation 0) is straightforward, and 
will be addressed elsewhere 
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